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A statistical field theory of particle production is presented using a gaussian func- 
tional in three dimensions. Identifying the field with the particle density fluctuation 
results in zero correlations of order three and higher, while the second order correla- 
tion function is of a Yukawa form. A detailed scheme for projecting the theoretical 
three-dimensional correlation onto data of three and fewer dimensions illustrates how 
theoretical predictions are tested against experimental moments in the different di- 
mensions. An example given in terms of NA35 parameters should be testable against 
future NA35 data. 
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1 Introduction 

Motivated by the search for a new exotic phase of matter predicted by QCD as well as a 
better understanding of nuclear matter at high energy densities, much experimental and 
theoretical effort is being directed towards studying heavy ions at high energies [|IJ. The 
fact that tens to hundreds of nucleons collide to produce several hundred particles makes 
it impossible to calculate analytically exclusive final states. Large-scale Monte Carlo simu- 
lations of microscopic processes and statistical methods have so far proven the only viable 
approaches to the problem. 

In the realm of statistical analyses, the measurement and theoretical treatment of one- 
particle observables such as cross sections, p_i_ distributions etc. have, after a period of dor- 
mancy, recently been supplemented by revived interest in correlations. Different theoretical 
attempts to describe correlations in various regions of phase space, classified under headings 
such as Bose-Einstein correlations 0, intermittency || |J and effects of collective flow ||, 
have met with some success but much remains unclear. Some recent experimental work sees 
the two as connected, saying that correlations of like-sign particles measured in intermittency 
analyses can be understood as a Bose-Einstein effect || , while conversely analyzing factorial 
moments in terms of the correlation integral and the four-momentum transfer Q 2 is equiv- 
alent to the Bose-Einstein analysis 0; in hadronic collisions, UAl has found a power-law 
dependence in their Bose-Einstein analysis ||. 

While the languages may be converging, this does not mean that there is as yet a solid 
theoretical basis for these correlations, especially those of higher order: it is still most diffi- 
cult to calculate measurable correlation functions from first principles. It is therefore helpful 
to explore simple effective field theories which ultimately should be derived from the under- 
lying theory of QCD. Thus the situation bears some resemblance to the pre-BCS era in the 
research of superconductivity. What we are looking for is a strong interaction analogy to 
the phenomenologically successful Ginzburg-Landau theory ||. 

The present paper attempts to understand correlations in heavy-ion collisions in terms of 
a statistical field theory based on a gaussian approximation of a suitably defined Ginzburg- 
Landau theory. This approach was first used by Scalapino and Sugar [[H]] who identified the 
field as a pion amplitude. In our case, the field is not an amplitude or particle density but 
a fluctuation from the mean density, an ansatz used previously |II[ |12| . It is motivated, as 
we shall see, by the fact that it matches experimental findings that in heavy-ion collisions 
only second order correlations are found while higher order correlations are negligible. In 
contrast to our previous work which was related to the rapidity variable only, all quantities 
are here calculated strictly in the three-dimensional phase space of rapidity y, azimuthal 
angle <fi and transverse momentum p± and then projected down for comparison with the 
data (see, however, the comment in Section VI concerning the proper choice of variables). 

The paper is organized as follows. In Section II, we remind the reader of the basic equa- 
tions for extracting true correlations and present the experimental evidence for neglecting 
correlations of order higher than two. In Section III, we recapitulate and develop the for- 
malism of statistical field theories in terms of a functional formulation and derive from a 
three-dimensional functional ansatz a general form of the second order correlation. We show 



in Section IV how one integrates out variables to find moments in one, two and three di- 
mensions and how this tests theoretical correlation functions. Using NA35 parameters and 
p± distributions, we show horizontal and vertical moments in all dimensions in Section V. 
A summary and discussion of some important issues conclude the paper in Section VI. 

2 Cumulants in heavy-ion collisions 

Traditionally, correlations were measured as a function of the distance between bins in phase 
space while keeping the bin sizes fixed. Following the proposal to look for a power-law struc- 
ture in the correlation function ||], a commonly used alternative has become the measurement 
of normalized factorial moments as a function of decreasing bin size while disregarding all 
distance information. 

The factorial moments F q (M) are constructed as follows. (For the purposes of this paper 
we shall stick to the coordinates rapidity, azimuthal angle and transverse momentum, but the 
formulation is true for all variables.) A given total interval Q tot = AY A0 AP is subdivided 
into M 3 bins of side lengths (AY/M, A0/M, AP/M). With n k i m being the number of 
particles in bin (k, /, m) and n™ = n(n — 1) . . . (n — q + 1), the "vertical" factorial moment is 

1 M (n [q] ) 



M3 HjfcLl ( n klm) q 

J_ ^ In klm U i d 3 x iPq {x 1 ...x q ) 



In klm d 3 xpi(x) 



The second equality illustrates how the factorial moment can be written in terms of integrals 
of the correlation function p q (flkim is the region of integration over bin k, I, m). Because for 
small bin sizes rikim becomes small and the relative error correspondingly large, an alternative 
definition is often preferred for three-dimensional analysis, the "horizontal" factorial moment, 

1 M (n [q] ) 

F h (M) = — V \ U klm/ 

= jkf3( ? -i) \p In klm Ujd x i p q (x 1 ...x q ) 



k,Lm=l 



In tot d3x Pi( x ) 



The latter form, while being much more stable, has the serious drawback that it is influenced 
by the shape of the one-particle distribution function p 1 . 

When one or more of the arguments of the correlation function become statistically 
independent it factorizes into lower order parts. To measure true particle correlations, known 
as cumulants, the trivial background must be subtracted. The first two cumulants are defined 
asUl 

C 2 (xi,x 2 ) = P2(x 1 ,x 2 ) - pi{xi)pi(x 2 ) , (3) 



C 3 (xi,X2,x 3 ) = p 3 (xi,x 2 ,x 3 ) - pi(x 1 )p 2 (x 2 ,x 3 ) 

- Pi(x 2 )p 2 (x 1 ,x 3 ) - pi(x 3 )p 2 (x 1 ,x 2 ) 

+ 2p 1 {x 1 )p 1 (x 2 )p 1 (x 3 ) . (4) 

By integrating these well-known relations over each bin, one can derive equations for inte- 
grated cumulants K v q = J Y\dxkC q / '(/ dxpi(x)) q in terms of the vertical factorial moments 
of Eq. (|Tp : the first two are 

Kl = F" - 1, Kl = F; - 3F" + 2 , (5) 



and so on for higher orders |T^]. Whenever there are no true correlations, these cumulants 
become zero. 

Analyses of the K q 's is now being carried out routinely in conjunction with the standard 
factorial moments as a function of bin size. Interestingly, it was found that in the case of 
heavy-ion collisions, there are only two-particle correlations: while K 2 is positive for all bin 
sizes, the values of K 3 , K4 and K§ have all been found to be consistent with zero. While 
initially this was found in terms of one- dimensional rapidity data, measurements by NA35 
in two and three dimensions have confirmed this fact ||L5|| . As an example we show in Fig. 1 
the third order cumulant derived from NA35 200 A GeV Oxygen-Gold measurements of the 



factorial moments [16]. (This data is given in terms of "flattened variables".) The higher 



orders K 4 and K 5 were also found to be consistent with zero. Corresponding findings for 



other nuclei and energies were published before |12], [T7|, [L8 



3 Statistical field theory 

For the large number of particles produced in high-energy heavy-ion collisions, a statistical 
theory of particle production is justified. A fruitful starting point has been an analogy first 
drawn by Feynman and Wilson fl9| , who pointed out that, by interpreting (y,4>,p±) as a 



"spatial" coordinate x, the final iV-particle phase space could be treated as if it were a gas, 
bounded by "walls" made up of the overall kinematic constraints. This analogy then carries 
over into a correspondence of the total cross section to a (grand canonical) partition function 
of a gas, while the n-particle cross section becomes the n-particle distribution function of 
the gas. It immediately leads to an identification of the probability of producing a secondary 
with a certain momentum (y,<f),p±) with the gas "density" pi at the corresponding point x. 
Botke, Scalapino and Sugar JL0[ ^(J utilized these ideas to write down a model for particle 
production. In analogy to ordinary statistical mechanics, where the density matrix governing 
the weights of states can be written in terms of the free energy, p oc exp(— (3(H — Q)), they 
defined a functional F[U] of a random field II, which governs the number of particles produced 
at point x via 

' (/VT (n 2 (£)> = 1 / VU e- F ^U 2 (x) , (6) 



a d 3 x ' ' Z 

where one identifies the particle density with II 2 and Z = JT>U exp(— .F[n]) plays the role 
of the partition function. Higher order correlation functions can easily be found by the 
appropriate functional derivatives. 



The essence of this ansatz is that all the physics is hidden in the form of the functional 
F[U] which reduces the greatly redundant information contained in the many microscopic 
degrees of freedom to just a few phenomenological parameters. (Botke, Scalapino and Sugar 
showed that their classical functional formalism is isomorphic to a system of quantized Boson 
fields with normal ordering [^(J . ) 

The analogy with the free energy permits the utilization of the Ginzburg-Landau expan- 
sion of the free energy near a second-order phase transition |T[| , which uses functionals of 
the type 



F[U(x)} 



d 3 x 



a(v x ri) 2 + /i 2 rr + A(ri 



T2\2 



(7) 



giving a minimum expectation value of the order parameter at zero above and nonzero below 
the phase transition. Applied to low-energy hadronic data by Scalapino and Sugar, it has 
recently been used for KNO scaling |22| and in an attempt to find a critical index for a 
phase transition to a quark-gluon plasma |23|] . Indeed, recent lattice gauge calculations 
suggest that, while a first-order transition is likely for the pure glue SU(3) gauge theory 
and for more than three quark flavors, there is increasing likelihood of a second-order phase 
transition for two quark flavors P5 |, and Wilczek [^S] has made first attempts at formulating 
such a transition in terms of a chiral order parameter. 

However, our aim here is more modest in that we do not specify the dependence of the 
constants a, fi and A entering the free energy functional in Eq. (|7|) on the dynamical param- 
eters of the heavy-ion reaction under consideration such as total center-of-mass energy, mass 
numbers and impact parameter. Instead, we treat them as phenomenological parameters to 
be determined from the experimental data and to be interpreted later on theoretically with 
the help of a more microscopic model as discussed in Section VI. Also we do not consider the 
definition of the functional of Eq. (f|) below an expansion valid for small values of the field $ 
only, and we cannot draw any conclusions whether the system undergoes a phase transition 
or not. An exception would arise when drastic and systematic variations in the phenomeno- 
logical parameters as deduced from one experiment compared to another are observed (see 
Section V). 

We define a random field $ as a function in the three-dimensional space spanned by 
(y,(j),p±). Throughout, p± will be implicitly divided by a constant scale V so that it is 
dimensionless. Since we are not looking for a phase transition, we omit the quartic term, 
A ~ 0, and start with the functional Jl2] 



r P f P/2 

F[$] = dy d 2 p ± 

JO J-P/2 



'f) 2 + « 2 (v^) 2 + /^ 2 



We note that F is not rotationally invariant but rather boost-invariant in the direction of the 
collision axis, i.e. tailored to the specific symmetry of the collision. The integration bound P 
is chosen the same for y and p± for simplicity but does not have to be the same in general. 
The expectation value of any function of $ is found by taking functional derivatives 
with respect to a source term, J ■ $, added into the integrand of F in Eq. (|8|) and thus 
into the exponential entering the functional integral for the partition function Z (take Xi = 



(Vi,(f>i,PAA)), 

1 5 k Z[F } J] 



Z 8J{x\) . . . 5J(xk) 



j=o 



(9) 



For the functional (H), we find the three-dimensional form of the two-point function 



(•<*)•&)> = ^^! , (10) 



where £ = a/p, and with p { = pj_j 



#- V(l/l -^2) 2 +Pl+P2-2piP 2 COs(0i-0 2 ) • (11) 

So far, we have not defined the field <J> in terms of physical observables. For reasons that 
will become apparent shortly, we define $(af) as the fluctuation at the point x of the particle 
density for a given event, pi(x), above/below the mean single particle distribution pi at that 
point: 

•(i) S £j*)-i. (12) 

Pi 0*0 
Experimentally, $ is the (normalized) difference between the event histogram and the event- 
averaged one-particle distribution for a given bin. 

This definition was previously used by Dremin and Nazirov [11| and in a one-dimensional 



context by Elze and Sarcevic ||12| . The rationale in the present and previous case is the same: 



all cumulants except the second order cumulant become exacly zero. To see how this comes 
about, we note that p q (xi, . . . , x q ) = (pi(xi) . . . pi(x q )), and with the definition ([12]) and the 
relations between the p g and cumulants C q one finds that the expectation values of $ can be 
written in terms of the reduced cumulants k q (xi, . . . , x q ) = C q (xi, . . . , x q )/px(xx) . . . Pi{x q ) 
as 

($(xi)$(x 2 )> = ^ 2 (xi,f 2 ), (13) 

($(fi)$(f 2 )$(f 3 )) = h(x l7 x 2 ,x 3 ) , (14) 

(${xi)&(x 2 )&(x 3 )Q(x A )) = A; 4 (xi,f 2 ,f3,f 4 ) (15) 

+ J] fc 2(^l,^ 2 )fc 2 (^3,^4) , 
(3) 

the sum running over 3 permutations in the arguments. On the other hand, the form of the 
functional (§) ensures that all expectation values of an odd number of $'s are zero, while 
for the fourth order, 

($(£i)$(£ 2 )$(30$(x 4 )> = £($(£i)$(£ 2 ))(<l>(£3)$0?4)>, (16) 

(3) 



so that &4 = also, and similarly for the higher even cumulants. These equations are actually 
just the relations between cumulants and central moments, see Eq. (3.38) in |T3|] . 



In summary, the specific form of the functional (|8]) and the definition of $ as a fluctuation 
field ensure that all cumulants of order 3 and higher are zero, while k 2 (x\, x 2 ) has the Yukawa 
form (|ig). 

We conclude this section with some comments: 

1. Of course with better statistics it may eventually become clear that there is some small 
residual cumulant of higher order; the presently large error bars would permit that. 
Were this to happen, the theory in its present form would have to be modified by, 
for example, the inclusion of interaction terms of higher order in <3>, with or without 
additional implications for the existence of a phase transition. 

2. The present model, where the higher order cumulants are exactly zero, is incompatible 



with the linking ansatz of Carruthers and Sarcevic |26[ used for hadronic collisions, in 
which higher order cumulants are products of k 2 according to k q = A q k 2 ~ , an ansatz 
which for certain values of the constants A q yields the negative binomial distribution. 



So far, there is no direct evidence for linking in heavy- ion collisions |17[; however, the 



present size of error bars does not permit any conclusions on this point. 

3. Thirdly, there has been some discussion whether cumulants of higher order can be 
identically zero consistently. Since this is a technical question with no consequences 
for our further development, we defer this point to the Appendix. 

4 Projecting down to lower dimensions 

Apart from the fact that all higher order cumulants are exactly zero by construction, the main 
result of our model is the Yukawa form for the second reduced cumulant, k 2 ex e~ R ^/R. This 
can be compared to data only after a suitable integration over its variables. In its present 
form, k 2 is not a power law of any of its variables. Setting cos 6 = 1, one can show analytically 
that J y d\yi — y 2 \ J P dpx dp 2 R~ 1 becomes linear in SySp, while numerical integration over 
all variables yields K 2 oc (5y5(p5p)~ a , with a = 1 within error, i.e. the integrated version 
resembles a power law. 

We now derive the detailed form of integrated cumulants K 2 for various dimensions in 
terms of the three-dimensional k 2 (xi,x 2 ). For this purpose, the "vertically" and "horizon- 
tally" normalized versions Eqs.([I|) and (0) have to be treated separately. 

The second order vertical cumulant is, in three dimensions, (always taking x = (y, 4>,p±)) 

1 M 

K v(6y,6^6p) = F 2 v -l = — £ mk,l,m), (17) 

k,l,m=l 

{n k i m (n k i m - 1)) - (rikim) 2 = ifw ^1 ^2 C 2 (x 1: x 2 ) 
( Uklm ^ 2 [Jn klm d3 x Pi(x)] 



with 

K*(k,l,m) 






Pi{xi)pi(x 2 ) 



^ kl d^ Pl {x)] 



In other words, the integration of k 2 to compare with data involves a correction due to 
the shape of the one-particle three-dimensional distribution function pi(x). Equation (|T8| ) 
as it stands is exact: a knowledge of the theoretical reduced cumulant k 2 can only be 
translated into a measurable factorial cumulant K 2 when the full three-dimensional one- 
particle distribution is known and taken into account (of course the same is true for r 2 = 
PilPxPi versus F 2 ). 

Knowledge of k 2 (x\,x 2 ) and pi(x) in three dimensions can also immediately be used to 
compare to factorial cumulant data of lower dimensions. For example, in (y, 0), the cumulant 
is K 2 (Sy, S(j)) = M~ 2 J2im ^(l, m) with the transverse momentum integrated over the whole 
window AP (cf. Section VI), 

^(i n (nim{ ni m - 1)) - (n lm ) 2 
K 2 {l,m) = 



(nirn) 2 

dyidy 2 /^ dfadfo j^dp 1 dp 2 - ^, (19) 

Ia m dy Iu t d<P J A p dp p x [x] 



n m -i^i mp 



and the cumulant for the rapidity only is K%(5y) = M" 1 X) m -^2( m )> with K^m) an in- 
tegral just like Eq. (|19"D but with J n being replaced by J A ^. Cumulants of other variable 
combinations are obtained analogously 

For the horizontal normalization, the three-dimensional cumulant K 2 = M~ 3 J2kim K 2 (k, I, m) 
consists of 

Ki(k, I, m) = ^m(n klm - I)) -(n klm ) 2 

where (N)a is by definition the number of particles within the experimentally defined total 
volume Qtot = AY Acf)AP . Care must be taken to define theoretical quantities such that 
they are normalized to this experimental domain. Keeping this in mind, we can write 



K\{k, I, m) = M 6 f d 3 x x d 3 x 2 k 2 {x 1 , x 2 )- 



pi(xi)pi(x 2 ) 



In tot d 3 x Pl (x) 

Projections onto two dimensions (y, 0) are then of the form of Eq. ( |19"D but with the pref actor 
M 4 and with Q t ot replacing the bin region integrals Qi, Q m , AP in the denominator. For 
one-dimensional data in y, the prefactor is M 2 and Qi becomes A0 in the numerator and 
Qtot i n the denominator integrals. Corresponding versions can be written down for K 2 (6(f)) 
and K^iSp). 

For horizontal factorial moments, it is important to remember that the simple relations 
between factorial moments and cumulants, Eq. (||), are not valid but that rather 

M 

F 2 h = K*+ £ (n klm ) 2 /(N) 2 n (22) 

k,Lm=l 



for three dimensions, with corresponding relations for lower dimensions. This has to be taken 
into account in an eventual comparison with data. 

With these relations it is thus possible, given any three-dimensional theoretical function 
&2 (or T2)i to compute factorial cumulants and moments for any combination of its variables. 
Doing this for different variables can serve as a strong test of the theoretical function as the 
moments probe its different regions. 

5 Comparison with the data 

So far, there are only two published sets of heavy-ion factorial moment data spanning several 
dimensions: one for 200 A GeV Sulfur-Gold collisions measured by the EMU01 collabora- 



tion |27|> |2l| and one- and two-dimensional data of the KLM collaboration |29fl. Both are 
unfortunately not suitable for our analysis, the former because of probable contamination 
of the (?7, 0) moment by gamma conversion |Kj , the latter because the data is binned with 



different M values in each variable. Both, being emulsion experiments, can measure only up 



to two dimensions. (An example of EMU01 analysis can be found in Ref. [31].) 

An analysis by NA35 of 200 A GeV Oxygen-Gold moments in all dimensions is in prepa- 
ration but not yet available; also, prospective Sulfur-Sulfur moments will provide further 
tests of the model |32| . In anticipation of such data, we have made a detailed study of mo- 
ments for NA35 experimental parameters and their measured O+Au transverse momentum 
distribution. Lacking the appropriate data, we have set the parameters to arbitrary but 
plausible values a = 2.0, £ = 1.0, hoping that exact fits can be made when experimental 
moments become available. 

Since all the equations of the previous section are exact, they can be used in their 
unabridged forms whenever the full three-dimensional pi(x) is known and enough CPU 
time is available. In making comparisons with data, it will however usually be necessary to 
approximate these exact forms, partly because three-dimensional one-particle distributions 
are not available, partly to save on computer time. We have therefore made the following 
approximations in our simulation: 

1. We factorize the three-dimensional one-particle distribution into its separate variables: 

p 1 (x) = (N) n g(y)h(<P)f(p A _), (23) 

where the factor (N)n ensures that the three distributions g, h and / are separately normal- 
ized over their respective total intervals AY, A$ and AP. (This factorization is known not 
to be true for some cases |£| .) 

2. The azimuthal distribution is taken as flat, h(<p) = 1/A$. 



3. We use the full experimental parametrization for f(p±) provided by NA35 |3"^ |; 

4. For the rapidity distribution, we use a gaussian parametrization with a = 1.32 |35 
and transform g(yi)g(y2) in the numerator with Y = (yi + y%)j2, y = yi — yi, to 



1 



g(Vi)g(V2) = —e- Y ^e-y'^ , (24) 



2trt 



and since k 2 (xi, x 2 ) is a function of \y\ — 2/2 j only, we have i?(AF) = J2 m In dY exp(— Y 2 /a 2 ) - 
j AY dY exp(— Y 2 /cr 2 ) as a constant prefactor, while the remaining gaussian in y is included 
in the numerical integration. 

5. As k 2 depends on the differences \y\ — y 2 \ and |0i — 02 1 only? we can transform to 
relative coordinates and integrate out the center-of-mass coordinate in these variables (the 
"strip approximation" |26f). 

With these approximations, we can derive the complete behavior of all cumulant moments 
both for the vertical and horizontal normalizations. In Fig. 2, we show the vertical cumulants 
for three-dimensional binning in (Sy,5(j),Sp), for the two-dimensional (5y,5(ft) and all three 
one-dimensional cumulants K 2 (Sy), K 2 (5(f)) and K 2 (5p±). The pseudo-power-law behavior 
is clearly visible for the three-dimensional case, while K 2 for two and one dimension show 



the familiar saturation caused by the projection process P5| . Fig. 3 shows the corresponding 
horizontal cumulants K 2 for the same set of variables. All curves of Figs. 2 and 3 are for 
a = 2.0, f = 1.0. 

A comparison of the two Figures shows that there is virtually no effect of the non- 
flat rapidity distribution, while the steep pj_ distribution has a large effect in raising both 
K%(5y,d<f),8p) and K%(8p) above their vertical counterparts. All cumulants must of course 
have the same value for M = 1, independent of their dimension, variable or normalization. 
This is also clearly illustrated in the Figures. 

Comparing to "real" experimental data, the procedure would be as follows: first, one 
takes the F 2 data of highest dimension and finds from these the corresponding K 2 from 
Eqs.(|5|) or fl22"D , depending on the normalization used. To this the free parameters a and £ 
are fitted using the appropriate formula of Section III. If the fit is disastrous, that of course 
is the end of the game. If it is reasonably good, further tests of the proposed cumulant 
k 2 (xi,x 2 ) are performed by plotting on experimental K 2 s of lower dimension and different 
variables the relevant theoretical formulae of Section III keeping a and £ fixed to their higher- 
dimensional best fit values. Of course, one can start out with lower-dimension best fits and 
find the corresponding curves for higher dimensions. This procedure is a nontrivial test of the 
theory because, for example, K 2 (5y) integrates over the entire regions A$ and AP and hence 
probes the long-distance behavior of the theoretical function k 2 , while the three-dimensional 
experimental K 2 tests only the short-range behavior of the theoretical k 2 . 

6 Conclusions 

To summarize our work, we studied integrated two-particle and higher order correlations 
in the form of factorial moments and cumulants and made predictions as to the general 
behavior of horizontal and vertical moments for the NA35 O+Au data. An exact comparison 
will depend on a fit of two parameters, after which everything else is fixed for all dimensions. 
It will be interesting to see how well our model will do with upcoming data from NA35 and 
hopefully other experiments. 

In particular, we applied a simple statistical field theory model with a gaussian "free en- 
ergy" functional motivated by analogy to the Ginzburg-Landau theory of superconductivity, 



10 



but without explicitly implying a (second order) phase transition. Our model, which is de- 



fined in terms of the field of density fluctuations, was first advocated in Ref. JLj to describe 
one-dimensional rapidity correlations as observed in high-energy heavy-ion collisions. Here, 
we extended the study to cover one-, two-, and three-dimensional correlations as measured 
through the intermittency analysis in the full phase space spanned by the variables rapidity, 
azimuthal angle, and transverse momentum of secondary particles. 

In our present formulation, the intrinsic scale V for the momentum cannot be deduced 
from the integrated dimensionless moments being compared with experiment; only direct 
measurements of the corresponding correlation functions, e.g. k q , would give enough detailed 
information for this purpose. Thus, also the mass parameter and correlation length of our 
model can only be determined as dimensionless constants. 

It should be stressed that (by construction) our model yields vanishing cumulant corre- 
lation functions of higher order by construction, fc g > 3 = 0, which agrees with all presently 
available heavy-ion data (cf. Section II). 

Working towards a more microscopic foundation of the Ginzburg-Landau type model, a 
comment about the right choice of kinematic variables seems appropriate: Since any attempt 
at deriving the effective three-dimensional statistical field theory from a more fundamental 
theory necessarily begins with a four- dimensional space-time formulation, i.e. leaving the 



Feynman- Wilson gas analogy |T9[ behind, three- or four-momenta conjugate to space-time 
coordinates are the natural variables. We therefore urge experimentalists to present their 
"intermittency" or correlation data in terms of three-momenta and with particle identifica- 
tion whenever this is feasible. Assuming azimuthal symmetry as before, the relevant variables 
for a more refined theoretical analysis seem to be p», p±, <fi. 

Of course, the difficult problem how to incorporate the effects of the very asymmetric 
initial conditions in p\\ and p± into the theoretical description remains as disturbing as before. 
Presently we attempted to include all available experimental information here by properly 
folding in at least the relevant one-particle distributions into the projection integrals for 
the lower-dimensional cumulant moments. An important consequence of this procedure is 
that, if our model fits a three-dimensional data set, then the projections onto two- and 
one-dimensional ones are predicted with no further freedom. 

It is also important to point out that functionals defined in different dimensions produce 
different results: K 2 calculated from a one- dimensional functional does not correspond to K 2 
obtained from deriving k 2 from a three-dimensional functional and then projecting down onto 



one dimension. Thus the formula obtained analytically in Ref. [I2| from the one-dimensional 
functional, K 2 = 2'j^ 2 [(Sy/^) — l + e~ Sy ^]/Sy 2 , does not hold in our present formulation (even 
though it may fit the data). 

Concerning the application of a statistical approach to energies higher than in current 
heavy- ion collisions, we have to deal with the expected increasing importance of hard partonic 
scattering events (minijets) |jj. They will help to populate more and more the high-p^ tails of 
the one-particle distributions, where a priori we cannot expect our model to apply. Whereas 
presently high momenta are effectively cut off by the p± distributions, this may become 
a problem at higher center-of-mass energies. Eventually a momentum cut-off has to be 
introduced to limit the analysis to the same region of soft physics implicitly investigated 
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here. A thorough treatment of surface effects in phase space, together with the asymmetric 
initial conditions [pTo[| , then becomes mandatory. 



Finally, we point out how a Ginzburg-Landau-type model may be based on a more 
fundamental field theory. Assuming that the observed multiparticle correlations and fluc- 
tuations arise in the dense hadronic phase of matter (after a possible hadronization phase 
transition), the first step of a derivation for this conservative scenario consists in choos- 
ing a phenomenologically satisfactory effective field theory of hadronic interactions such as 
the sigma model. (This presents, of course, an alternative to the view that the relevant 
fluctuations stem from partonic showering processes which might seem more natural for 
high-energy e + e~-collisions.) Keeping the high excitation energy of the compressed mat- 
ter in mind, one tentatively begins with a field theory at finite temperature, neglecting for 
simplicity the difficult non-equilibrium aspects of heavy-ion collisions. Aiming at a three- 
dimensional Ginzburg-Landau statistical picture one is led to integrate out the (imaginary) 
time-dependent modes from the path-integral representation of the (sigma model) field the- 
ory. This has so far been achieved in a one-loop approximation for an arbitrary scalar field 
theory yielding a temperature-dependent three-dimensional effective action^. It represents 
the infrared limit of the originally chosen theory in accordance with ideas on dimensional 
reduction and can be used as a starting point to calculate the correlation functions studied 
in the present paper (see Section III). 

Clearly, much more can be done here to include non-perturbative effects if one wants 
to investigate, for example, effects of a phase transition possibly contained in the chosen 
field theory. Formally integrating out the fields of the three-dimensional effective action 
in favor of the fields squared ("densities"), one obtains a Ginzburg-Landau "free energy" 
functional with its parameters, cf. Eq. (0), given as temperature-dependent functions of the 
renormalized coupling constants of the original field theory. 

Having outlined these further steps in the development of our model, we conclude that the 
study of correlations and fluctuations in dense hadronic matter seems a promising approach 
to further our understanding of strong interactions. In particular, the investigation of one- 
particle observables alone, which are nicely reproduced by a wide selection of event generators 
PJ, in principle cannot provide a sufficiently detailed knowledge of the most interesting (and 
least understood) soft physics aspects of QCD. 
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Appendix: 

Finding multiplicities from cumulants 

When all factorial cumulants are known, it is usually possible to derive from them the 
multiplicity distribution P n . This is done via the factorial moment generating function 

Q(X) = £(1 - \) n P n , (25) 

n 

which can also be expanded in terms of the unnormalized cumulants f q = JC q as 



Q(A) = exp 



E 



■A)« 



g! 



./; 



(26) 



For our model, only the first two cumulants are nonzero, f\ = n, the average total multiplic- 
ity, and f 2 = K 2 n 2 . With appropriate transformations and the identity P n = (l/n\)(—d/d\) n Q(\)\\ = 
one can derive the multiplicity distribution in terms of Hermite polynomials 



Pn = -xhm nl ^- h+hl2 H n 



TV. 



*(/ 2 - /i)A/2/ 2 



(27) 



from which, for example, Pi = (/i — f 2 )Po- Clearly Pi becomes negative for f 2 > ft, 
rendering the multiplicity distribution invalid. In terms of normalized factorial moments, 
this requires K 2 < n^ 1 for Eq. fl2~T|) to be valid. All odd P„'s are similarly dependent on the 
sign of the factor (/i — / 2 ), while the even ones are positive. For most heavy-ion data, the 
average multiplicity is very large, meaning that the multiplicity distribution fl2"T|) cannot be 
trusted for such cases. 

Clearly, something is amiss. The resolution of this dilemma is found on returning to 
the identification of the normalized cumulant functions k q with the expectation values of the 
fluctuation fields, Eqs.(fL3|) if- In this identification, the difference between factorial cumulants 
and ordinary cumulants is neglected, which for the large multiplicities of heavy ion collisions 
is negligible. Mathematically, however, there is a difference, and this shows up in the above 
dilemma for the derivation of P n from the factorial cumulants. 

If we work with ordinary cumulants, however, we can unambiguously derive the multi- 
plicity distribution as follows: the characteristic function 



cj.it) = £ e' 4 "P n 



:28) 



ra=0 



is also the sum of ordinary cumulants n q (e.g. K\ = (n) = n; n 2 = (n 2 ) — (n)' 



(p(t) = exp 



Eyltj K q 
q =i y- 



which for our truncated cumulant set is 

(f)(t) = exp itn — t 2 K 2 /2 



(29) 



(30) 
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and since (up to a normalization constant) the multiplicity distribution is 

/oo 
e~ itn <P(t) dt , (31) 

-oo 

we get 

P n oc exp — (n — n) 2 /2/€ 2 , (32) 

just an ordinary gaussian multiplicity distribution which is defined perfectly well. 



Figure captions: 



Figure 1: Third order cumulant K3 as a function of the number of bins M for NA35 OAu 
data in (y, 4>,p±) ||16| . Cumulants of higher order are also compatible with zero. This fact is 
confirmed in analyses in terms of other variables and different colliding nuclei. 



Figure 2: One function ^2(^1, ^2) determines all: Theoretical vertical cumulant moments K\ 
for various dimensions, for fixed parameters a = 2.0, £ = 1.0, incorporating the experimental 
NA35 rapidity and p± distributions for 200 A GeV O+Au. 



Figure 3: Theoretical horizontal cumulant moments K\ for the same fixed parameters 
and NA35 distributions as in Fig. 3. The effect of the p± distribution is clearly visible. 
Comparison with NA35 data requires conversion to horizontal factorial moments F% and a 
fit of a and £. 
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